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This  technical  note  provides  the  forecaster  the  necessary  information  to 
understand  the  physics  used  in  these  products.  Understanding  the  strengths  and 
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INTRODUCTION 


The  replacement  of  the  GEOVOR  computer  at  the  German  Military  Geophysical  Office 
(GMGO)  near  the  end  of  1986  made  possible  the  use  of  better  atmospheric  models 
due  to  increased  computing  capacity.^ 


Specifically,  the  synoptic  scale  models  7LPE-CM  (hemispheric,  gndpoint  spacing 
381  km)  and  7LPE-LFM  (North  Atlantic/European  window,  gridpoint  spacing  127  km) 
were  replaced  by  a  more  powerful  model.  These  models,  used  by  the  GMGO  since 
1982,  were  limited  by  the  model  physics  and  numerical  methods  available  in  1980. 
Advancements  in  these  fields  gave  modelers  the  power  to  upgrade  these  models.''-. 


"  Due  to  the  existing  cooperation  between  the  DWD  and  the  GMGO  in  the  area  of 
model  development,  an  agreement  was  reached  to  implement  the  DWD  hemispheric 
model  BKF  at  the  GMGO  in  a  version  specifically  adapted  to  the  GEOVOR  computer 
system.  This  new  version  is  called  the  BKFG  (G  for  GMGO),' 


After  extensive  program  adaptation  and  testing,  the  BKFG  was  available  for 
routine  use  on  GEOVOR  and  replaced  the  7LPE  as  the  standard  GMGO  model  on  3  May 
1988. 

The  BKFG  is  a  9-layer  primitive  equation  model  with  an  effective  horizontal 
resolution  of  180  km.  The  computer  code  corresponds  to  the  DWD  version  of  the 
BKF  as  of  March  1986.  The  quality  of  the  forecast  is  highly  dependent  on  the 
initialization  fields.  Deviations  from  DWD-products  can  therefore  result  due  to 
different  analysis  procedures,  the  use  of  climatology  instead  of  analysis  in  the 
surface  model,  and  a  different  data  supply  due  to  a  different  ‘cutoff'  time.  The 
analysis  procedures  are  continuously  being  improved,  and  the  climatology  fields 
are  replaced  by  analyses  in  the  working  process. 

In  spite  of  these  limitations,  a  verification  test  of  the  BKFG  forecasts 
performed  by  the  GMGO  during  a  6-month  period  revealed  significant  improvement 
in  the  important  products  over  the  7LPE-LFM  model.  Regardless  of  the  somewhat 
poorer  horizontal  resolution,  better  vertical  resolution,  more  detailed  model 
physics,  and  the  avoidance  of  window  model  boundary  value  problems  account  for 
the  improved  products.  For  a  detailed  description  of  the  BKFG  model,  the  reader 
is  referred  to  the  GMGO  Int  .nal  Report  Nr.  79137. 
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BKFG  COARSE  MESH  ANALYSIS  &  FORECAST 
NUMBER  OF  GRIDPOINTS:  79  X  79 
GRIDPOINT  SPACING:  254  km 


BKFG  FINE  MESH  FORECAST 
NUMBER  OF  GRIDPOINTS:  65  X  65 
GRIDPOINT  SPACING:  127  km 
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BKFG  BASIC  EQUATIONS 


a.  PROGNOSTIC  EQUATIONS: 

1.  Horizontal  equations  of  motion: 

f-"1 


-$*•»(*>) 


First  Law  of  Thermodynamics: 
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♦  =  H  *D(8) 
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Moisture  Conservation  Equation: 


Tendency  Equation: 

Mr-  .  .  RT0 

61  ■UoP» 


DIAGNOSTIC  EQUATIONS: 

1.  Hydrostatic  equation: 


^=-0)1'Ke|,.°6om 


2.  Continuity  equation: 

-  m’(pi,pL 
£>P  \6x  6y 

TERMS  EXPLAINED: 


Vv,  V2  Covariant  horizontal  velocity  components 

Vertical  velocity  in  the  X.Y.P  coordinate  system 
6  Potential  temperature 

q  Specific  humidity 

H  Sources  and  sinks  for  heat 

R  Sources  and  sinks  for  humidity 

D  (Y>  Turbulent  flux  of  variable  Y  due  to  internal  friction  and  boundary 
layer  phenomena 

m  Scaling  factor  (stereographic  plane) 
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BKFG  MODEL  PHYSICS 


SURFACE 

OROGRAPHY 

BOUMDARY  LAYER 

COrVECTlOR 

CLOUDS 

PRECIPITATIOM 

RADI  ATI  OK 


Composition:  land,  ocean,  ice,  and  variable  snow  cover.  Two 
layers  with  prognostic  variables  temperature  and  soil  water 
content . 

Vertical  walls  (mountains)  inhibit  flow.  Air  is  lifted  over 
the  barrier. 

Combination  of  Prandtl  and  Ekman  layers  allowing  calculation  of 
boundary  layer  fluxes  using  surface  roughness  and  Richardson 
number  as  the  determining  parameters. 

Parameterization  of  deep  convection  by  the  parcel  method. 
Vertical  velocity  is  calculated  from  the  lifting  force  and 
determines  the  intensity  of  the  convective  exchange.  The 
following  are  considered: 

-  slow  sinking  of  the  surrounding  air 

-  entrainment 

Parameterization  through  linear  dependence  on  relative  humidity 

Consideration  of  stratiform  and  convective  precipitation. 
Condensate  is  not  stored  in  clouds  ar.d  immediately  falls 
through  lower  model  layers  (where  it  can  evaporate) ,  reaching 
the  surface  as  rain  or  snow  depending  on  the  temperature. 

In  the  troposphere  (950-300  mb) : 

-  heating  by  shortwave  radiation  depends  on 
--  sun  angle 

--  water  vapor  concentration 

-  cooling  by  long-wave  radiation  (1.8  +  0.017  *  T  ®C/day) 

--  reduced  below  clouds 

--  increased  above  inversions 

Upper  surface  layer: 

-  shortwave  heating  depending  on  sun  angle,  cloud  cover,  albedo 

-  long-wave  cooling  depends  on  cloud  cover  (counter-radiation) 
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BKFG  HUMEBICAL  METHODS 
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Octagonal,  79  x  79  grid  points  with  254km  grid 
spacing  (at  60  N)  using  stereographic  projection. 

Muller  grid  ('staggered'):  all  variables  except 
defined  at  the  grid  point;  2  types  of  grid  points 
(red/black)  diagonally  shifted  by  1/2  grid  point 
with  distribution  of  variables  depending  on  time 
step;  -*  effective  grid  spacing  160  km. 

x,  y,  p  coordinate  system  (0  <_  p  <_  lOOOmb) 

9  layers  for  Vi,  v2  ,  6,  q 
2  surface  layers 

Explicit,  centered  (leapfrog  scheme) ,  At  =  1/14 
hour.  Averaging  over  2  starting  periods  (to  damp 
out  high-frequency  initial  waves)  and  over  double 
time  steps  (to  avoid  diverging  solutions) 
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c . 

Bottom  (p  *  1000) 

w  =  0  (kinematic  boundary  condition) 
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BKFG  INITIALIZATION  FIELDS 
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BKFG  CALCULATED  PBODUCTS 
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BKFG  WE AX MESSES 


1 .  General  weaknesses 


-  The  BKFG  is  limited  by  the  octagonal  grid  requiring  significant  data 
smoothing  south  of  30°N. 

-  Grid  point  models  cause  truncation  errors  due  to  the  data  only  being  input 
and  calculated  at  a  specific  point. 


2.  Model  physi cs 


-  Homogeneous  surface  composition  doesn’t  take  into  account  the  difference 
between  urban  vs  wooded  or  rural  areas.  Temperature  and  humidity  errors  can 
occur . 

-  Missing  analysis  of  ice  coverage  can  cause  problems  for  the  radiation 
scheme  when  calculating  sensible  heat  flux  and  albedo. 

-  Uncertainties  in  boundary  layer  parameterization  due  to  missing  Ekman 
spiral.  Turning  of  the  winds  isn’t  considered. 

-  Missing  cloud  water  phase  in  precipitation  scheme  can  cause  over  forecast 
of  precipitation  since  the  model  precipitates  out  any  water  as  soon  as  it 
reaches  the  condensation  level. 

-  Highly  simplified  radiation  scheme  for  surface  and  free  atmosphere 
(950-300  mb)  is  used  due  to  the  amount  of  time  it  would  take  to  run  a  more 
complex  scheme.  Products  would  be  delayed  a  few  hours. 


Numerical  Methods 


-  Explicit  time  integration  requires  small  time-step  making  computation  time 
longer . 

-  Leapfrog  scheme  yields  phase  errors  with  short  waves  causing  them  to  be 
too  slow. 

-  Poor  vertical  resolution  in  the  boundary  layer  and  in  the  vicinity  of  the 
tropopause . 

-  Artificial  lateral  boundary  conditions  require  significant  smoothing  (no 
conservation  of  mass) . 

-  Lower  boundary  condition  represented  in  the  x,y,p  coordinate  system  not 
the  most  efficient  but  it  reduces  interpolation  error. 

4.  Initialization 

-  No  differentiation  during  time  averaging  between  fast  synoptic  waves  and 
very  long  'noise’  waves  with  the  same  frequency. 

-  Octagonal  model  area  does  not  permit  'non-linear  normal  mode’ 
initialization. 


VEBIFICATIO!  COMPABISON 
BKFQ/7LPE 


The  verification  procedure  teste  forecast  versus  observed  changes  with  respect 
to  the  average  values  of  the  parameters  in  question.  It  is  valid  for  the 
Central  European/East  Atlantic  window  and  gives  reliability  averages  for  the 
testing  period  of  1  Nov  87  -  29  Feb  88.  Statistics  show  that  the  BKFG  did 
significantly  better  than  the  LFM  through  the  72-hour  point.  Notice  that  both 
models  tend  to  get  worse  near  the  72-hour  panel.  Care  should  be  taken  when 
using  these  extended  forecasts. 


I.  Correlation  Coefficient  R  (Phase  error) 
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II.  Root  Mean  Square  (RMS)  Error  Q  (Amplitude  Error) 
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BKF/BKFG  VERIFICATION 


1, 


500  mb 

CORRELATION  COEFFICIENT 


BKF/BKFG  VERIFICATIOB 


TRANSMITTED  PRODUCTS  WITH  DESCRIPTIONS 


Following  are  the  BKFG  fields  routinely  transmitted  over  EURDIGS  to  all  2WW 
units.  Descriptions  of  the  parameters  found  on  each  field  follow  the  panel. 


Surface  pressure  in  millibars  and  850  mb  temperatures  in  C.  H  represents  high 
pressure  while  T  represents  low  pressure.  CAN  depict  areas  of  cold  and  warm 
air . 


13 


—  )RELTGP  500/1000  12Z  +  12HBKFG  24.05.89 


( 


500/1000  mb  thickness.  Heights  in  geopotential  meters. 
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(—  )  GEOPOT 
)  TEMP 


700  HPfl 
,700  HPR 


1 2Z  +  12HBKFG  24.05.89 
1 2Z+  12HBKFG  24.05.89 


700  mb  geopotential  heights  in  decameters  and  ’00  mb  temperature  in  °C. 


( —  ) GEOPOT  500  HPfl  12Z  +  12HBKFG  24.05.89 
(--  ) TEMP  500  HPR  12Z  +  12HBKFG  24.05.89 


500  mb  geopotential  heights  in  decameters  and  500  mb  temperature  in  °C . 


( —  JRELTOP  500/850  00Z+  QOHBKFG  22.06-89 

22.06.89  1 :30.0  MIO 

EDZX  NR  K  .56 1 


I 


)REL  FE  700  HPfl  12Z  +  12HBKFG  24.05.89 


700  mb  relative  humidity  in  X.  */-  refer  to  areas  of  maximum  and  minimum 
relative  humidity. 


TECHNICAL  NOTE  DISTBIBUTION 


AWSTL ,  SCOTT  AFB  IL  62225-5008  5 

AWS/DOZ ,  SCOTT  AFB  IL  62225-5008  1 

AFGWC/DO ,  OFFUTT  AFB,  NE  68113-5000  1 

3350  TECHTG/TTGU-W,  STOP  62,  CHANUTE  AFB  IL  61868-5000  1 

1  WW/DN,  HICKAM  AFB  HI  96853-5000  1 

2  WW/ 

OL-B/OIC,  2  WW,  APO  NY  09305  1 

DET  6.  2  WW/CC.  APO  NY  09126  1 

OL-A,  DET  6,  2  WW,  APO  NY  09050-5000  1 

DET  13,  2  WW/CC,  APO  NY  09122  1 

WL,  APO  NY  09094-5000  1 

7  WS/DON,  APO  NY  09403-5000  26 

28  WS/DON,  APO  NY  09127-5000  8 

31  WS/DON,  APO  NY  09136-5000  17 

3  WW/DN,  OFFUTT  AFB  NE  68113-5000  1 

4  WW/DN,  PETERSON  AFB  CO  80914-5000  1 

5  WW/DN,  LANGLEY  AFB  VA  23665-5000  1 

7  WW/DN,  SCOTT  AFB  IL  62225-5008  1 


